Spontaneous symmetry breaking of gap solitons in double-well traps 
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We introduce a two-dimensional model for the Bose- Einstein condensate with both attractive and 
repulsive nonlinearities. We assume a combination of a double-well potential in one direction, and 
an optical-lattice along the perpendicular coordinate. We look for dual-core solitons in this model, 
focusing on their symmetry-breaking bifurcations. The analysis employs a variational approxima- 
tion, which is verified by numerical results. The bifurcation which transforms antisymmetric gap 
solitons into asymmetric ones is of supercritical type in the case of repulsion; in the attraction model, 
increase of the optical latttice strength leads to a gradual transition from subcritical bifurcation (for 
symmetric solitons) to a supercritical one. 
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I. INTRODUCTION 

The Gross-Pitaevskii equation (GPE) provides a pow- 
erful model for studying the mean- field dynamics of Bose- 
Einstein condensates (BECs) Important examples 
are the prediction of ID gap solitons (GSs) in a self- 
repulsive condensate trapped in a periodic optical-lattice 
(OL) potential 2]. This was realized experimentally in 
an ultracold gas of *^Rb atoms confined in a cigar-shaped 
trap H, and the prediction of the Josephson effect in a 
BEG [j.It was subsequently observed in a condensate 
trapped in a macroscopic double- well potential [B]- In 
contrast to hitherto realized Josephson systems in super- 
conductors and superfluids, interactions between tunnel- 
ing particles play a crucial role in a bosonic junction. The 
effective nonlinearity induced by the interactions gives 
rise to new effects in the tunneling. In particular, an- 
harmonic Josephson oscillations were predicted d, S Hi j 
provided that the initial population imbalance in the two 
potential wells falls below a critical value This dy- 

namic regime can be well explained by means of a simple 
model derived from the GPE, which amounts to a sys- 
tem of equations for the inter-well phase difference and 
population imbalance. The nonlinearity specific to the 
BEG also gives rise to a self-trapping effect in the form 
of a self-maintained population imbalance. 

One-dimensional dynamics of a BEC in potentials com- 
posed of two rectangular potential wells were studied in 
several papers [11]. Stationary states with different pop- 
ulations in the two wells are generated by symmetry- 
breaking bifurcations from symmetric and antisymmetric 
states, for attractive and repulsive nonlinearity, respec- 
tively [i,iTo[. A natural 2D extension of the double-well 
configuration is a dual- channel one, with the potential 



featuring the two wells in the direction of x, which are 
extended into parallel troughs along the y axis [T^ . [l3| . 
In the case of an attractive nonlinearity, this setting 
may naturally give rise to dual-core solitons, which are 
self-trapped in the y direction (similar to the ordinary 
matter- wave solitons created in a single-core trap .14]), 
and are supported by a double-well structure in the per- 
pendicular direction. Furthermore, if the nonlinearity is 
strong enough, or else the tunnel coupling between the 
troughs is weak, the obvious symmetric dual-core soliton 
may bifurcate into an asymmetric one. This was demon- 
strated both in the full 2D model [l^], and in its ID 
counterpart, which replaces the 2D equation by a pair 
of one-dimensional GPEs with coordinate y, while the 
tunneling in the x direction is ap pro ximated by a linear 
coupling between the equations [ij]. In fact, the lat- 
ter model resembles the standard one widely accepted in 
nonlinear optics to describe dual-core nonlinear optical 
fibers and asymmetric solitons (isl . [l6j . In a similar way, 
the double-well potential may be uniformly extended in 
two transverse directions, giving rise to a 3D structure 
based on a pair of parallel "pancakes" . 

If the dual-channel potential in 2D geometry is com- 
bined with an axial optical lattice, which runs along both 
potential troughs, it is natural to consider a dual-core gap 
soliton in the self-repulsive BEG filling this structure. 
In Ref. [l^, this was done using the above-mentioned 
approximation which replaced the corresponding two- 
dimensional GPE by a pair of linearly-coupled ID equa- 
tions. It was demonstrated that a symmetric gap soli- 
tons may be stable in this case, and never bifurcate, 
while asymmetric solitons are generated by a symmetry- 
breaking bifurcation from antisymmetric ones. Similar 
results (including the emergence of asymmetric gap soli- 
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tons carrying intrinsic vorticity) where obtained in the 
2D extension of the model. This model pertains to 
the above-mentioned "dual-pancake" structure In 
nonlinear optics, asymmetric gap solitons were studied 
in models of dual-core fiber Bragg gratings, which also 
amount to systems of linearly coupled ID equations p^ . 

The prediction of symmetry breaking for matter-wave 
solitons in a setting combining the transverse double- well 
potential and a longitudinal optical lattice in experimen- 
tally relevant conditions makes it necessary to study the 
full 2D model (especially for the stability of the emerg- 
ing asymmetric solitons) for both repulsive and attrac- 
tive condensates, which is the purpose of the present 
work. Parameter regions admitting asymmetric solitons 
will be predicted by means of the variational approxima- 
tion (VA) |16] . These results will be verified by numerics. 
The character of the symmetry-breaking bifurcations for 
the dual-core solitons will also be identified (we obtain 
a gradual transition from a subcritical bifurcation to a 
supercritical one with increase of the OL strength). 

The paper is organized as follows. The model and the 
VA are introduced in Sec. II. In Sec. Ill we analyze the 
symmetry-breaking bifurcations in both attraction and 
repulsion models, and Sec. V concludes the paper. 



II. THE MODEL AND VARIATIONAL 
APPROXIMATION 

The normalized form of the GPE for the mean-field 
wave functions 5* in 2D geometry is 

i-^t = -(1/2) (*.. + ^yy)+[U{x) + a\^f + pcoa (2y)] 

(1) 

where a — +1 and —1 for the self-repulsive and self- 
attractive BEG, and pcos(2y) represents the longitudi- 
nal optical lattice potential. The transverse double-well 
structure is taken as 

fO, < L/2 and |x| >L/2 + D, 

'^^'')-\ -Uo, L/2 < \x\ < L/2 + D, 

with D, Uq and L being, respectively, the width and 
depth of each well, and the width of the barrier between 
them, see Fig. [1] below. 

Stationary solutions to Eq. ([1]) are assumed in the form 
^{x,y,t) = e~*'^*$(a;, y), where the real function ^{x,y) 
satisfies the equation 

/i$ + (1/2) ($,, + ^yy) - U(x)^ - + p cos(2y)$ = 0. 

(3) 

It can be derived from the Lagrangian, 

istat = y y dxdy - (1/2) ($2 + $2) _ 

-U{x)<P^ - {a/2)<i>* + pcos{2y)<i>^] . (4) 

To apply the VA, we follow Ref. [12] and adopt an ansatz 
consisting of two distinct parts. First, inside each po- 
tential trough, i.e., at \x^ {L D) /2| < D/2, the trial 




FIG. 1: (Color online) The shape of the quasi-one-dimensional 
double-well potential, U{x,y). The wiggles indicate quasi-lD 
lattice along y. 



function is 



^ / ^ A , xt{L + D)/2 
^± \x,y) = A± cos I TT ) exp 



D 



2W^ 



(5) 

where A± and W are three variational parameters. This 
expression implies different amplitudes and a common 
longitudinal width, VF, of the wave-function patterns in 
both troughs. In the x direction, the ansatz ([5]) emulates 
the ground-state wave function in an infinitely deep po- 
tential box, which vanishes at the edges of the trough, 
see Fig. [TJ In the y direction, the ansatz approximates 
the self-trapped soliton by a Gaussian profile. Outside 
the troughs (at |a;| > L/2 + D and \x\ < L/2), the ansatz 
also follows the pattern of quantum mechanics, in the 
form of a superposition of exponential wave functions: 



^{x, y) = ^A± exp ( 
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(6) 

with the same amplitudes A± and width W as in Eq. ([5]). 
The ansatz is not continuous at the edges of the troughs; 
however, comparison with numerical findings (see Fig. [2] 
below) clearly suggest that the VA can be used despite 
this local discrepancy. 

Substitution of expressions P and dll) into Eq. g]) and 
integration produce the following simplified Lagrangian, 
in which contributions from the exponentially decaying 
functions in the outer region, \x\ > L/2 + D, are ne- 
glected, the contribution from the optical lattice poten- 
tial is taken into account only inside the troughs, and 
the Thomas-Fermi approximation in the x direction is 
adopted, i.e., term — (1/2)<I'2 in the Lagrangian density 
is omitted: 
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We now define N± = (3/4\/2) A^W , and 



A EE (2/1?) v/^2^exp -y^(L + i:) 



N+ + N_ N+- 



4^A ' 



(8) 



(9) 



(10) 



The numbers of atoms trapped in the two troughs are 
proportional to the respective partial norms of the wave 
function, 
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dy 



±(D+L/2) 



±L/2 



dxmx,y)y 



2V2^ 



-DN± , 



(11) 

hence defined in Eq. (jlOp . measures the population 
imbalance. In this notation, the Lagrangian ^ simplifies 
to 



-Lcs 



(12) 



^eN N VX N^ + iy^ r— .J -w- 

2 8W^ 2 W 2'^ 

(13) 

with s = +1 and —1 for the configurations of the anti- 
symmetric and symmetric types (with < and 
A+A- > 0, respectively). 

Our Lagrangian gives rise to variational equations 
dL/dW = dL/di' = dL/dN = 0: 



N + 2(tV\ {N^ + v'^)W - ApNW^e-^ = 0, (14) 

) = 0' (15) 



W 



iV2 



2VXN 
W 



2sXN 



pe 



e. (16) 



Equation has two solutions: 1^ — 0, which corre- 
sponds to symmetric or antisymmetric solitons, and 
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(17) 



for asymmetric ones. Comparison of typical asymmetric 
and symmetric solitons, found from a numerical solution 
of Eq. ([3]), with their counterparts predicted by the VA, 
is presented in Fig. [21 

For symmetric and antisymmetric solitons, Eqs. 
and (|16p , with v — 0, are tantamount to equations that 
were derived, by means of the VA, for solitons in ID 
models with a periodic sinusoidal potential and attractive 
or repulsive nonlinearity [iqI . [20| . In particular, in the 
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FIG. 2: (Color online) The top and bottom panels demon- 
strate examples of cross-section profiles, along y = 0, of sta- 
ble asymmetric and symmetric gap solitons in the model with 
repulsion, as obtained from a numerical solution to Eq. (|3]) 
and predicted by the variational approximation (dashed and 
continuous lines, respectively). Parameters of the double- well 
potential are L = D = 1, Uo = —0.7 (repulsive case) and 
p = 1. Norms of the asymmetric and symmetric solitons are, 
respectively. A'' = 0.52 and 0.34. The asymmetry parameter 
for the former soliton, see Eqs. (|10|) . is v — 0.34. 
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FIG. 3: Three dimensional version of the situation pictured 
in Fig. [2] . The density scale is represented by the intensity 
of the print. Horizontal stipes represent double well structure 
and modulation illustrate the lattice 



latter case (for <t = +1) a known fact is that solutions 
exist only for p > p'"' = e^/16 ~ 0.462 (in fact, this 
constraint predicts, with high accuracy, the edge of the 
first finite bandgap in the linear spectrum induced by the 
OL [131). Results for asymmetric solitons are presented 
in the next section. 



III. ASYMMETRIC SOLUTIONS 



A. Equations for the bifurcation point 



According to Eq. (|15[) . asymmetric solutions exist in 
two cases: a ^ s — +1 (repulsion, with the asymmet- 
ric branch bifurcating from the antisymmetric one), or 
a = s = —1 (attraction, with the bifurcation from the 
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symmetric branch). Elimination of in Eqs. ((T4)) and 
by means of Eq. (|17p yields a system of equations 
for and VF: 



(18) 



1 2^/\N 2sVXN 



W 



W 



pe 



Taking into account definitions (O and (jlOp , solutions to 
Eqs. (|18p depend on parameters L, Z?, Uq, and p. 

At the bifurcation point, v — 0, Eq. P7|) yields 
N = VXW, hence Eqs. ^ generate a system of two 
equations for two coordinates of the bifurcation point, fj, 
[via relations HO]) and ([9])] and W: 



l + 2aXW^ = ipW^e'^ 



(19) 



+ 2 (cr + s) A - pe" 



Without the OL, i.e., for p = (the case considered in 
Ref. [12]), the first equation in (fT9|l gives the bifurcation 
point at TV = l/-\/2. To obtain explicit results in the 
model with p ^ 0, one can start with an obvious solution 
to Eqs. at X = = N = 0, p ^ p(o) (j.g^j^ii ^(0) ^ 



e2/16), U = U^^' = 1/16, and W = W^"'' = V2. This 
solution, which has N — is, by itself, trivial, but a 
nontrivial one can be obtained as an expansion around 
it. 



B. The model with self-attraction 

Consider the attraction model corresponding to cr = 
s = — 1. Then, straightforward analysis of Eqs. (fT9|) for 
small Sp = p—p^^^ and SUq = Uq—U^^ demonstrates that 
the bifurcation of symmetric solitons (which pertain to 
s = — 1, see above) may occur at two values of the norm. 



N = 



1 



2\/2 



{e-Hp + 5Uo) 



2 - - {e-'^6p + SUo) ± y^l5e-^6p- 5Uo 



(20) 



the respective value of the width being W « 
\/2 [l - {e-^Sp + SUo) /4] ■ Note that the second term 
in the square brackets in Eq. (|20p is a small correction 
to 2, the main correction given by the last term, which 
demonstrates that theoretically there may be two differ- 
ent bifurcation points. Obviously, expressions (|20p are 
meaningful, i.e., the bifurcation takes place, if 



(eVl5) SUo <Sp< --e^SUo 



(21) 



(in other words, SUq must be negative, while Sp may have 
either sign). Numerical calculations imply that only the 
lower value of N is valid. 




FIG. 4: (Color onIine)A set of numerically found bifurcation 
diagrams in the model with attraction, showing degree of a 
asymmetry of dual-core soliton, i^, as a function of the soli- 
ton's total norm, A'^, see Eqs. pop . The diagrams pertain to 
fixed values of parameters of the transverse double- well con- 
figuration, L — D — 1, Uo = —0.7 (attractive case), while the 
strength of the longitudinal optical-lattice potential gradually 
increases. One can check from the analysis of Eq. (|19p that 
the turning points are at = and ±A''/\/3. One can clearly 
see that the supercritical bifurcation will turn into subcritical 
bifurcation with increase of the optical latttice strength. 



A set of bifurcation diagrams in the attraction model, 
in the form of v{N), i.e., curves showing the asymme- 
try of the dual-core solitons versus the total norm, was 
generated by a numerical solution of the full system of 
Eqs. (HH]). The set is displayed in Fig. 21 where a 
noteworthy feature is the transition from the subcritical 
shape (backward-directed one), which is a characteristic 
of the attraction model without the longitudinal OL [13[ 
(as well as to the model of dual-core optical fibers [ISj). 
to the simpler supercritical (forward-directed) shape at 
sufhciently large values of OL strength p. Note that 
the symmetry-breaking bifurcations of dual-core solitons, 
studied in systems of linearly-coupled GPEs including the 
attractive nonlinearity and OL potential [H, [13] , as well 
as in the system of linearly-coupled fiber Bragg gratings 
[isl ]. are of supercritical type too. The physical signifi- 
cance of the subcritical bifurcation is that it allows bista- 
bility of the solitons (the coexistence of stable symmetric 
and asymmetric ones) in a limited interval of values of 
N. 



C. The model with self-repulsive nonlinearity 



In the case of the self- repulsion, i.e., a — s — +1, the 
expansion of Eqs. ^9^) predicts the following values of 
the norm at which asymmetric gap solitons may bifur- 
cate from the antisymmetric ones (recall antisymmetric 
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U=-4, 1=1, d=l 



■ p-0.6 
. p=0.8 



FIG. 5: (Color online) A set of bifurcation diagrams for gap 
solitons in the model with repulsive nonlinearity, for L = D = 
1, Uo = 4, and a set of different values of the OL strength, p. 



solitons corresponds to s = +1) 



2 - - {e-^Sp + SUo) ± v/l5e-2<5p - 6Uq 



(22) 



pliance with results obtained for the models based on 
linearly coupled GPEs with the optical lattice potential 
and repulsive nonlinearity [13, [l3| ■ These bifurcation di- 
agrams exist only for p > = 6^/16, because, as said 
above, at smaller values of the optical lattice strength the 
VA does not predict antisymmetric GSs that might give 
rise to a bifurcation. 



IV. CONCLUSIONS 

We have introduced a 2D model for self-attractive and 
self-repulsive BECs, which combines a double-well po- 
tential in the transverse direction, and a periodic poten- 
tial along the longitudinal coordinate. The analysis in- 
volved symmetry-breaking bifurcations for dual-core soli- 
tons. Systematic results were obtained by means of the 
variational approximation, which was verified by numer- 
ical results. In the case of a repulsive nonlinearity, the 
bifurcation is of supercritical type, while in the model 
with attraction an increase of the optical lattice strength 
leads to a gradual transition from subcritical bifurcation 
to a supercritical one. This is an important result. 



where the notation is the same as in Eq. ([20]) for the 
attractive model. This expression predicts the bifurcation 
in the following region [cf. Eq. ((2T|) in the attraction 



model] : 



Sp < SUo < 15e Sp, which implies dp > 0, 



while SUq may be both positive and negative, in contrast 
with the case of the attraction model, that demanded 
SUq < 0, while allowing Sp to take either sign. Once 
again numerical calculations imply that only the lower 
value of N is valid. 

A typical set of bifurcation diagrams in the repulsive 
model is displayed in Fig. [51 It is seen that the bifur- 
cation generating asymmetric gap solitons from the an- 
tisymmetric ones is always of supercritical type, in com- 
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